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Abstract 

We have used Monte Carlo simulations to observe the magnetic behaviour of Ising thin-films with cubic lattice 
structures as a function of temperature and thickness especially in the critical region. The fourth order Binder 
cumulant is used to extract critical temperatures, and an extension of finite size scaling theory for reduced geometry 
is derived to calculate the critical exponents. Magnetisation and magnetic susceptibility per spin in each layer 
1-^3 ' are also investigated. In addition, mean-field calculations are also performed for comparison. We find that the 
Q ■ magnetic behaviour changes from two dimensional to three dimensional character with increasing thickness of 
the film. The crossover of the critical temperature from a two dimensional to a bulk value is also observed with 
both the Monte Carlo simulations and the mean-field analysis. Nevertheless, the simulations have shown that the 
critical exponents only vary a little from their two dimensional values. In particular, the results for films with up 
to eight layers provide a strong indication of two dimensional universality. 
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I. INTRODUCTION 



Studies of the dimensional crossover of magnetic properties from two dimensional {2d) to three dimen- 
sional {3d) character in magnetic multi-layers have currently attracted much interest. An understanding of 
the physical properties of solids as the dimensionality is reduced is of both technological and fundamental 
importance In particular, studies of ultrathin magnetic films have revealed a number of novel phe- 

nomena that would not have been expected in either 2d or 3d systems. Of particular interest is the critical 
behaviour of magnetic thin-films, for which the dimensionality d is not well established. It is interesting 
to consider how magnetic properties such as the magnetisation m, magnetic susceptibility x ^^nd critical 
temperature Tq are affected by the finiteness of the system in the direction perpendicular to the film. 

From both theoretical and experimental studies, the critical temperature Tc is found to increase with 
film thickness. Nevertheless, how the thickness / of the film relates to the dimensional crossover of other 
magnetic critical properties from 2d to bulk 3d is not fully understood. Most studies imply that magnetic 
films belong to a 2d universality class regardless of the thickness. This is commonly understood to be a 
consequence of the finiteness of the films along the out-of-plane direction. Close to the phase transition 
point (critical temperature Tq) the correlation length ^ is constrained by thickness and allowed to expand 
only in the in-plane direction. This is certainly a 2d-\ike phenomenon. However, well-known experimental 
studies of thin-films of nickel jsj have shown contrary evidence where a dimensional crossover of the critical 
exponent jS from 2d to 3d has been found. Hence, to clarify this discrepancy, we consider the use of Ising 
thin-films to study such a phenomenon in cubic lattices, that is simple cubic (sc), body centred cubic (bcc) 
and face centred cubic (fee). 

The Ising model represents a ferromagnet with infinite uniaxial anisotropy. It is useful because strong 
magnetic anisotropies are common in ultrathin ferromagnetic films and monolayers. In addition, it turns 
out that 2(i anisotropic Heisenberg systems become Ising-like near Tq. This was shown by Binder and 
Hohenberg using Monte Carlo simulations and has been proven rigorously by Bander and Mills |5i|using 
renormalisation group analysis. Also, from the experimental point of view, experiments on nickel j3[ and 
iron thin-films p have shown 2d results in agreement with Ising 2d critical exponents. These all suggest 
that the Ising model is a handy tool to study magnetic thin-films. 

Previous studies of Ising thin- films have been by means of the high temperature expansion 0, 0|, 
renormalisation group j9|, variational cumulant expansion (VCE) [i3 and Monte Carlo Il^ methods. 
However, most of them have concentrated on the study of the critical temperature and shift exponent in 
the sc lattice. Few have investigated the effective critical exponents in order to determine how they vary as 
a function of temperature away from Tc- As far as we are aware, nothing exists in the literature reporting 
the investigation of critical exponents as a function of thickness at the critical point. This is mainly due to 
a presumption that the films are of the 2d class as well as difficulties in the analysis due to finite size effects. 
Hence, in this study, we aim to give a more complete picture of the magnetic phase transition in thin-films 
especially at the critical point. We have investigated how the magnetic properties vary as a function of 
temperature and thickness by means of Monte Carlo simulations and mean-field theory. We have also 
calculated the critical temperature and investigated its convergence to the bulk limit. Most crucially, 
we have developed an extension to finite size scaling analysis to observe how the critical exponents vary 
as a function of thickness. In outline, we firstly describe the numerical calculations and the methods 
we have used. Secondly, we show the evolution of the magnetic properties as a function of temperature 
and thickness. Then, we report the critical properties of thin-films in terms of critical temperatures and 
exponents. Finally, we discuss our results and compare the characteristic critical exponents with those 
found in experiments. 
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II. METHODOLOGY 



The starting point for the study is the nearest neighbor Ising hamiltonian 

H = -JY.S,S,, (1) 

<ij> 

where the spin Si takes on the values ±1 and the sum includes only first nearest-neighbor (Inn) pairs. 
Helical and free boundary conditions are used in our simulations using this hamiltonian for the in-plane and 
out-of-plane directions respectively. We use units of J/ks and J for temperature and energies respectively 
with the magnetisation per spin defined as m = Si where is total number of spins. 

The simulations are carried out for sc, fee and bcc films of size N = L x L x I where L x L represents 
the number of sites in each layer of the film and / is the number of layers in the film i.e. its thickness. 
We vary L from 64 to 128 (in steps of 8) with / ranging from a monolayer (bilayer for bcc films) to 20 
layers. The spin configurations of the films are updated using the highly efficient Wolff algorithm 
to minimise the effect from statistical errors ari sing from correlation time 0,0]. The random number 
generator (drand48) is chosen carefully 0, 17, [l8| and we ensure that different seed numbers have no 



significant effect on the results. For each simulation, roughly 3, 000 x spins are updated before the 
system is deemed to have reached equilibrium. From this point, 5 x 10^ independent configurations are 
used to calculate the expectation of magnetisation per spin < m >= -r^ X^i'""" l^^iL where tmax is total 
number of measured configurations, and the magnetic susceptibility x = > — < \m\ >^) ( 

f3 = J/ksT). In a similar fashion, the layer-dependence of these magnetic properties, mk and Xk-, where 
/c is a layer index, are calculated in order to observe the surface effects upon the magnetic properties. For 
the investigation of the critical behaviour of the films, we locate the critical temperatures (Tc) via the 
fourth order cumulant Ul [3] 

3 <m"^ >^ 

At T = Tc, Ul should be independent of L, i.e. for differing sizes L,L' {Ul' /Ul)t=Tc = 1- Owing to finite 
size effects, we need to plot Tc{b) (L = bL') against (ln6)~\ and extrapolate the results to the infinite limit 
i.e. (In 6)^^ — > To maximise the efficiency of this Tc calculation, for each thickness, we perform a 



single long simulation at a temperature Tq and use the histogram method [20|, |21| to extrapolate Ul to 
a temperature nearby in order to find the cumulant crossing point. This temperature is chosen from the 
temperature at the peak of the magnetic susceptibility curve for the L = 128 system, and then around 1 
to 4 million spin configurations are used to create the histograms. To exclude those data obtained from 
temperatures too far from the simulated temperature Tq, the range of extrapolation \T — Tq\ is restricted 
by the criterion \U{T) — U{Tq) \ < (Te [25] where U =<E>, the average of the energy, and a^; is a standard 
deviation of E at Tq. 

The next step is to extract the critical exponents from our finite size results. To do this we have 
developed a finite size scaling method appropriate for films in which l/L ^ 1. The purpose is to find how 
magnetic propert ies m and y scale with the size L and / of the systems. The basic finite size scaling ansatz 
(e.g. (UHll in HIS) rests on an assumption that only a single correlation length ^ is needed to 
describe the critical properties of thin-films. As detailed in Appendix ^ this gives the following general 
form for how the magnetisation and susceptibility scale with L 

<m{T,l)> = L-^/"m(L^/'^t,0, 

X{TJ) = L^^Z-^yxiL'/^l), (3) 



where {•y/i^Y = •y /u + 2 — d. 7, /3 and z/ are the critical exponents associated with x, ^ and ^ respectively, 
X and m are scaling functions for a given / and t = 1 — ^ is the reduced temperature. These scaling 
functions for a range of L should collapse onto a single curve with the correct critical temperature and 
exponents. The exponent 1/u can be extracted from the derivative of the cumulant with respect to L at Tc 
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owing to its variation with system size as L^/" [l^ • Since we have assumed that there is a single ^ oc a 
single V is to be expected. Note that if equations Q correctly encapsulate the nature of magnetic critical 
behaviour in films, we can extract the exponents (3/v and (7/1^)' by plotting logm and logx against logL 
at Tc- Before this can be done with confidence it is necessary to demonstrate the validity of the assumption 
of a single correlation length and thus the form of equations This can be done by establishing the 
following :- 

1. According to eqn. 0, at Tc, if we plot logx {'^n) against logL we should get a straight line. 

2. We should find that the sum of (7/z/)' + 2(/5/z/) is always equal to 2. This again can be derived using 
the assumption of a single correlation length and '-y /v + 2l3 /u = d ji^ i.e. 

f^V + 2^ = l + 2-d + 2^ = 2, (4) 

3. The clearest demonstration comes from a direct examination of the scaling functions themselves, 
eqn. Q. With correct exponents in the critical region, a scaling function for any L but a particular 
/ should collapse onto a single curve. This will also confirm the reliability of the critical exponents 
extracted from our simulations 



III. RESULTS AND DISCUSSION 

A. Magnetisation and Magnetic Susceptibility 

The magnetisation m and susceptibility x profiles for various film thicknesses I and system sizes L are 
obtained. The crossover of magnetic behaviour from 2(i-like for the monolayer (bilayer in bcc) to 3(i-like 
for films with 20 or more layers is found. For example. Fig. ^,b shows the magnetic phase transition 
for a system of 128 x 128 x I spins. The transition point moves from 2c? to values with increasing 
film thickness. The layer resolution of these magnetic properties is also investigated. The influence of the 
surface on the magnetic properties can also be seen. For instance. Fig. \^ and d show the layer-dependence 
of mfc and Xk for the 128 x 128 x 10 system. The surface layer has the lowest magnetisation magnitude, 
whilst the magnetisation per site increases towards the interior of the films. Similarly, for x, representing 
the fluctuation of m, the strongest magnitude is for the innermost layer and the weakest is at the surface. 
The fc*^ and (Z — fc)*'' layers have the same magnetic properties owing to the spatial symmetry and the 
isotropy of J. 

The layer variation of the magnetic properties can be explained in terms of finiteness of the thin-films 
along the z direction together with the Monte Carlo updating algorithm. In a simulation run, once a seed 
spin is chosen in order to select a group of spins to be updated with the Wolff algorithm, it will have 
a greater chance of aligning those neighboring spins which are located towards the interior of the films 
rather than those towards the surface. Consider the probability for a site getting updated: 

^ M 

D{t,j,k) = -Y.s'{z,j,k), (5) 

^'^ n=l 

where s'{i,j, k) is the number of times the site at location k) is in a Wolff update, and M is the total 
number of clusters forming (flipping) in a simulation. Our simulations show that, in the ferromagnetic 
phase, D at the surface layers always has the lowest magnitude, except, at very low T where the chance of 
a site getting updated is very high and the same anywhere. Thus, this implies that once a cluster of spins 
is formed, the group has its spins preferentially located in the inner layers rather than near the surface 
layers. When the group is flipped, the inner layer will always have the high magnetisation magnitude, 
but at the surface layers, there will be more un-flipped spins behaving as a buffer which results in smaller 
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magnetisation magnitude. Similarly, for the susceptibility, the buffer in the surface layers helps to reduce 
the magnetic fluctuation leading to a smaller susceptibility than for the inner layers. 

We find the same qualitative trends by using a mean-field approximation (see appendix The 
evolution of the magnetic properties from 2d- to 3(i-like is again found e.g. for sc films in Fig. |2ti^,b. 
The layer-dependence of magnetic properties is also evident. As expected, the magnetisation per site is 
smallest on the surface layer and highest on the innermost as can be seen, for example, in Fig. |2fc,d which 
show results for a 10-layered Ising sc film. 



B. Critical Temperatures and Shift Exponents 

We calculate the critical temperatures of the films from our Monte Carlo simulations using the cumulant 
method. The results are summarised in Table E] We find a change from 2d to bulk values as / is increased. 
Fig. |2t shows evidence of such a dimensional crossover for the Ising thin-films. Both our 2d and bulk results 
agree very well with the exact Ising 2d results and previous Ising 3(i studies [21. 30j . As the number of 
layers is increased, the critical temperature moves towards the 3(i value owing to the increase of the average 

1 r 

exchange interaction energy. On similar grounds, for Z > 4, we find T^^{1) < T^^^{1) < Tlp^{l). The 
mean-field theory calculations show the same qualitative trends (Fig. Eb) but, as expected, the estimated 
Tp's are consistently higher than those from the Monte Carlo work. Moreover an analytic expression for 
Tc as a function of film thickness has been made on the basis of a mean-field Tc examination [s^l • This 
is given as, 

^ .. + 2..cos(./(i+l)) ^ 

Zq + 2zi 

where zq and zi are number of nearest neighbors in the same plane and its adjacent plane respectively. 
We find our mean-field results to comply with this expression. A summary of results for Tc from both 
Monte Carlo and mean-field calculations is presented in Table H] and Fig. El 

It is revealing to examine the evolution of the thin-films' critical temperatures from the monolayer to 
bulk "id limit in terms of a power law 

Tc(oo) 

Here Tc{oo) is the bulk critical temperature and the shift exponent A has a value between 1.0 and 2.0 
depending on the spin model used and the type of calculation. For thick films A is expected to be 
|27| . A better fit for films of a range of thicknesses / is given by jl^ and has the form 



1 1 



Tc{l) Tc(oo) 



1 + 



where /q, /' and A' are all adjustable parameters. Similarly, A' should tend to as / oo. 

We use eqn. (jS)) to fit the Tc(/)'s arising from our Monte Carlo calculations. If it turns out that our 
T(7(/)'s are accurate and that the fit is a useful one, we should find Tc{oo) in agreement with the T^^^ 
we calculate separately. Results of the fit are shown in Table |nj and there is less than a 1 percentage 
difference between Tc{oo) and T^^^^ for all three types of films. As can be seen in the Table ITTl however, 
even for / = 20, A' is not close to the expected large I value, To elucidate further the evolution from 

2d- to 3(i-like behaviour we rearrange the power law of eqn. ((7j) and define 



Ae.(/.) = - log /Mt^), (9) 




and tabulate Aefr(/i) with /j for both our Monte Carlo simulations (Zj G {4, 6, 8, 10, 15, 20}) and mean-field 
calculations (/« ranging from 4 to 20 layers) in Table HH The Aeff 's should converge to asymptotic shift 
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exponents Xjv^'^ when / tends to infinity. A linear least squares fit between Aefr(/i) and enables us 
to obtain Aefr(cxD) which is also given in the Table. (The mean-field values of Aefr(cxD) for all structures 
all have values of 2 since the mean-field v is well-known to be 1/2.) For the Monte Carlo simulations' 
results, it is gratifying to find Acfr(oo) to be close to Xjv^'^ which we obtain from our separate bulk Monte 
Carlo simulations. This gives good support to the contention of universality and the asymptotic behaviour 
contained in eqn. ()7l8p . 

C. Critical Exponents 

Referring to equations Q, we notice that the critical exponents can be extracted from plots of logm 
and logx, as well as log ^j^, against logL at Tc, provided that linear relations are found. This is actually 
the first test of the validity of our single correlation length, ^, assumption. Results from our simulations, 
at T(7, indeed show very good linear relationships between logm, logx and log with logL for all our 

thin-films' thicknesses and structures. We note here that oc Li^l^ . An example is shown in Fig. 0] 
which presents very good linear fits for 10-layered sc films. Now, the exponents (7/1/)', and Xjv can 
be extracted from the slopes of the linear least square fits. They are listed in Table HI For the second test, 
we perform the summation of the computed exponents + 2^. We have found that, in each case, this 
sum has a value of 2 within error bars as shown in Fig. El According to eqn. (^, this passes the second 
test for the assumption of a single ^ being sufficient to describe the critical behaviour of the thin-films. 
For the last test, we consider the scaling functions in eqn. Q. We calculate these scaling functions for 
the 10-layered sc films with L = 10,20,40,60,80,100,150 and 200 using the exponents extracted from 
L = 64 to 128. The results are shown in FiglHti^jb. For large enough L (L = 60 to 200), excellent data 
collapses occur. This confirms that our critical exponents are perfectly reliable. However, in the same 
figures, for small L (around L = 10 and 20), data collapsing is not found. This is what we must expect 
since eqn. dH)) originated from eqn. ()A6|A7|1 (in appendix EJ under the condition Z -C L. A collapse of 
scaling functions should not be observed if L approaches /. For comparison, we also present the scaling 
functions for these 10-layered sc films using the 2d monolayer critical exponents in Fig|ni:;,d. We can see 
that the data collapses are not quite as good. This also confirms the reliability of our results. 

Since our results pass all of the three tests above, it can be understood that the single ^ assumption 
is valid and eqn. Q is very useful for extracting critical exponents from film systems at Tc- However, it 
is unfortunate that we cannot extract the 7/z/ out of (7/z^)' from our results. This means that we cannot 
investigate the dimension ci as a function of thickness. Nevertheless, we find that the exponents (7/z^)' 
and jS/u, for / = 1 up to / = 20 films, are quite close to their 2d values as shown in Table H] and Fig. El 
In particular, for thin- films (/ < 8), the exponent values seem to be almost identical with the 2d results. 
On the other hand, for thicker films (/ > 10), a weak dependence of the exponents on / is found. This is 
somewhat reasonable since, in the critical region, although the correlation length ^ is 2(i-like in the sense 
that it can expand freely only in the xy plane, it is still a function of thickness. Suppose that at a reduced 
temperature t close to zero, any film has an in-plane correlation length ^y. Then, the total magnitude of 

the correlation length ^ ^ \/ (^||)^ + different thickness / will also be different. The thicker the film, 
the more different from 2d and the closer to 3d it will be. However, thick films are beyond the scope of 
this study and we find that, for thin-films with thickness around 8 layers and thinner, the exponents are 
very close to the 2d values. Thus, we feel that it is quite safe to conclude that thin-films belong to the 
2d class. In addition, our results show that, at a good level of agreement, the critical exponents from all 
structures represent the same universality class. 

Finally, we compare our results for the exponent f3 with those available from an experiment with Ni 
films 0. As shown in Fig. [71 the experimental results are close to ours only for films thinner than 4 
layers. This is not surprising because the anisotropic Heisenberg magnet becomes Ising-like only for very 
thin-films. However, instead of a sharp dimensional crossover of (3 at around 5 to 8 layers 12], we found a 
trend which suggests that a dimensional crossover from 2d- to Sd-Ising should occur for much thicker films. 
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This is perhaps indicated by the shght increase oi j3/v towards its bulk value from / = 10 onwards. This 
issue may also be partly answered in that the dimensional crossover of /? in may really be a transition 
from 2(i-Ising to 3(i-Heisenberg instead of from 2c?- to 3(i-Ising in our study. So, a direct comparison may 
not be allowed. Nevertheless, although Ising and Heisenberg films may be quantitatively different, they 
should qualitatively share same characteristics since, in the critical region, the divergence of ^ in both 
models is 2(i-like. So, it is strange that /? in ||3| should change its value so abruptly with very few layers. 
To find more answers, we consider 2| (Schilbe et al,1996) which claimed 'the author of [3] neglected the 



dependence of the (effective) critical exponents away from the critical point on the reduced temperature 
in their evaluation of the experimental data'. In other words, the range of temperatures 10"^ < t < 10'^ 
used for the power law fit in Q may not belong to the asymptotic behaviour (critical region), and the fit 
may lead to somewhat dubious results. So, to investigate this closely, we follow who defined 

_ dlogm 
ologt 

and study how this varies with t in thin-films. 

From Monte Carlo simulations in sc Ising thin-films for temperatures below but close to Tc, results of 
logm against logt for L x L x 10 sc films are shown in Fig. |H| The figure shows how, close to Tc, for each 
L, the correlation length ^ notices the finite size L, that is where the finite size effect becomes important. 
Roughly, depending on the value of u (see Table E)), ^ oc t^'^ will realise its finite limit at a temperature 
around logt = — ^ log L. We have also plotted Pcs against logt for 100 x 100 x I films in Fig. IHl It is found 
that the pattern of the jSes curves is in agreement with the rough investigations in |1^. Since L = 100 
is used in our calculation, the PeS are reliable (with minimum finite size effect) only down to logt ~ —2. 
However, this range of logt does cover most of the range —3 < logt < —1 used in the power law fit in 
and shows a quite interesting phenomenon. It appears that this range of temperature is not in the critical 
region because f3ea for each film is not constant, but peaks at a certain temperature. Outside the critical 
region, the correlation length grows as we increase the temperature towards Tc- As long as ^ is smaller 
than the film thickness, the film tries to behave as the 3d bulk system and (3 grows somewhat towards the 
bulk value /^^^ ^ 0.3258 |2l|. However, at some temperature ^ becomes comparable with the film thickness 
and is only allowed to expand in the in-plane direction, that is 2d behaviour. Thus, f3 changes its trend 
and the resulting decrease results in a peak. We can also notice that the thicker the film, the closer the 
peak can grow towards the bulk value. Of course, if one tries to perform the power law fit in this range of 
temperature, a sudden change of PeS will be observed. 

We may conclude that the range —3 < logt < —1 used for the power law fit in is outside the critical 
region and the power law fit should not work well. On the other hand, it must be emphasised that our 
exponents are extracted at the critical temperature by means of finite size scaling. Our exponents are the 
real critical exponents. This explains why the behaviour of /3 in |3i] should not relate to our results. 



IV. CONCLUSION 

We have studied the magnetic behaviour of Ising thin-films in sc, bcc, and fee structures using extensive 
Monte Carlo simulations. We have found the dimensional crossover of both the magnetisation m and the 
magnetic susceptibility x from 2d- to 3(i-like with increasing film thickness. The layer-components of m 
and X have also been presented. The surface layers are found to have the lowest magnitudes whilst the 
innermost layers have the highest, a trend which is derived from the free boundary at the surface. (This 
behaviour is also found in our mean- field study.) We track how the films' critical temperatures, Tc, evolve 
from 2d to 3d values with increasing film thickness. Results for square monolayers agree well with the 
exact Ising 2d Tc and our bulk 3d results are also in good agreement with earlier studies in the literature. 
We are also able to fit our calculated Tc's for films of varying thickness, /, involving 3 parameters (shift 
exponents). This expression could be examined in the limit of infinitely thick films and the limit agrees 
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very well with that calculated directly from 3d simulations. This fitting expression could also be used on 
the mean-field calculations. 

We examine the critical regime of these systems in detail and extract critical exponents. For this 
purpose we develop a finite size scaling method for films. It is based on one assumption that a single 
correlation length ^ = ^{i) is sufficient to describe the critical behaviour of Ising films. The validity of this 
assumption is successfully verified. From our results we find a very weak variation of the critical exponents 
with respect to /, the film thickness. For thin-films (/ < 8) the exponents are essentially the same as 2d and 
from this it can be implied that thin-films fall into the 2d class. For thicker films, however, (/ > 10) a weak 
/ dependence is noticeable because / is thick enough for the correlation length to distinguish the geometry 
of the films from that of a simple 2d lattice. As an example we present figures of the scaling functions for 
10-layered sc films which are constructed by using both their own 10-layer critical exponents and also the 
well-known 2d exponents. For large enough L (layer extent) we find good collapsing (L-independence) for 
the thickness-specific exponents which underscores the accuracy of our results. This collapse degenerates 
for L I which confirms that our finite size scaling analysis is valid only for L ^ I. When we use the 2d 
critical exponents in the scaling functions a slightly poorer data collapse is found which points to the real 
disparity between thin-films and 2d exponents. 

In comparison with the experimental results described in our results for the thinner films bear up 
well. A direct comparison, however, is not possible. This is because the experimental data from nickel 
films on a tungsten substrate can be interpreted to show a transition from 2d-lsmg to 3(i-Heisenberg rather 
than from 2(i-Ising to 3(i-Ising which is the only possibility for the model we have studied here. There is 
also the issue that the experimental data used to make a power law fit are taken from temperatures that 
are outside the critical regime. 
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APPENDIX A: FINITE SIZE SCALING 

We assume that the films' extent shows up when the correlation length growing along the z direction 
{^z) is about the same size as /, but C,x and ^y, the in-plane x and y correlation lengths, are still smaller 
than L and approach L when T — > Tc- So, for a film in the critical region, we propose a hypothesis that 
a single ^ is enough to describe the critical behaviour i.e. C,x = C,y = ^ with .^^ = Z as a constant. This 
implies ^ = ^(/). From this we modify a technique described by Binder and Wang 33] . The aim is to find 
out how thermodynamical properties such as m and x scale with L in a finite size lattice in the critical 
region where the film thickness / is a constant. 

Consider the susceptibility 

Lx Lz Lx Ly 

ksTx = T T T X] 51 H H <m(xi,?/i,zi)m(a;2,l/2,22)> - <m>^ (Al) 

X y z xi=l yi=l zi=l X2=l y2=^ Z2=l 

where we set the lattice spacing as 1. We use periodic boundary conditions along the x and y directions. 
Changing the index Z2 to -|- 2; under a condition that Lx and Ly 3> 1, we can write 

l rLx rLy L^-zi 

ksTx^ — / dx / V V <m{0,0, zi)m{x,y, zi + z)> - <m>'^ . (A2) 
Jo Jo ,^^1 

Next, sup pose that the two-point correlation function for an isotropically shaped system at Tc takes the 
form (3^ 

Gf(r,Tc) = <m(0)m(r)> - <m>2|^^ ~ -J^ ~ (^^^ + y^ + z^p''-^+'^^/^ (A3) 
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(where r = |ri — r2|, ci is the dimension of the system, and rj is a critical exponent) and that this form is 
also valid for anisotropically shaped systems like films, we rewrite the susceptibility 



ksTcxiTc) - Idx' Tdy'Y. E + ^ + , (A4) 

where we have scaled the variables x L^x' and y Lyy'. In the thin-film structure, we set = Ly = L 
and Lz = I, so that 

kBTcxiTc) - / dx' dy'Y. E h^" + ?/" + 7i • (A5) 

As 2 C [!-/,/-!] and if we choose / < L, ^ < 1, this yields Eii=i Ei=l-2i(a;'^ + 1/'^ + I^) ~ {x'^ + y''^) ■ 
So the susceptibility becomes 

kBTcxiTc) - KL'Y-'^ fdx' ['dy'{x'' + y")-^'-'^^^/' 

JO JO 

oc V'/'^+^-'^l oc (A6) 



where the scaling relation 7/1^ = 2— 77 (e.g. [2^) has been used. Note that for a particular thin-film system, 
/ is a constant so the only variable in the function is L. Consequently, in thin-films, if / -C L, x scales with 
L in the same way as isotropic-shaped systems, but with a different exponent (7///)' = 7/// + 2 — rf. For 
< m >Tci 'we may assume that it is of the same order as the root mean square magnetisation < >Tc^- 
So, 

1/2 



2 1/2 

<m>Tc oc <m >r^^ 



ksTcXiTc 



c) 



m 

oc L-^/^ (A7) 



where again we have used 7/(2z/) — d/2 = —(3/1/ j2^. Note that since ^ = ^(/) in our hypothesis is a 
function of thickness /, every exponent calculated in this way will also be a function of thickness. 

APPENDIX B: MEAN-FIELD APPROXIMATION 

In a mean- field study of thin- films, the average Weiss field is layer-dependent. We write down the 
expression for the free energy for an Z-layered films and obtain the equations of state from its minimisation 
with respect to the layer-resolved magnetisations {k is a layer index). The /-coupled equations are 
solved numerically to extract the equilibrium magnetisations m^. The magnetisation per site in the k*^ 
layer is given as rrik = -Pfc.t ~ Pk,i where Pk^tl) probability of a site being occupied by an up (down) 

spin. Including the interaction with an external field {—hj^i Si) the interaction energy U is the sum over 
all layers (fc = 1, ■■■,/) of terms Ei^ which are written as 

Ek = {ZoJkkml + ZiJfe,fc+imfcmfc+i(l - 5k,i) + Jfc,fc~i"^fcmfc_i(l - 4,i) + 2hmk] , (Bl) 

where Zq and Zi are the number of nearest neighbors to a lattice site in the same layer and one of its 
adjacent layers, and N\\ is the number of sites in a layer. The free boundary condition necessitates the 
presence of the factors 1 — 6k^i and 1 — 6k,i in the equation. Similarly the entropy S can be written as 
S = EL=i Sk where 



Sk = -ksNii [Pk,r lnPfc,t + (1 - Pfc.r) ln(l - Pk^)] 



-ksNw 



l + TTifc l + nik l-nik. l-ntk 
m 1 m 



(B2) 
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Thus, by minimising the free energy F = U—TS with respect to nik i.e. = 0, we obtain the equihbrium 
magnetisation in layer k by solving the following /-coupled equations for /^layered films 

- (ZoJi,imi + Zi Ji + h) + — — In = 0, 

2 LI — mi. 

-{ZoJk,kmk + ^1 Jfc,fc+i^fc+i + ^1 Jfc,fc-i^fc--i + h) + In ] =0, = 2, 1 

2 LI — mfc. 

- (Zo J;,/m„ + ZiJ,^/_imi_i + /i) + In | =0. (B3) 



Tc can also be extracted from eqn. ()B3j) for zero external field. Close to Tc, the {nii} are small and 
In ~ 2mj and Tc can be extracted by finding the temperature for which the set of equations Am = 0, 
where A is a / x / matrix with elements 

Aj = {ksT - ZoJij)6ij - ZiJ.j {(1 - 6i^i)6ij-i + (1 - 6i^i)6ij+i} , (B4) 

and m is a / X 1 column matrix,{mi, ■ ■ ■ , mi}, is satisfied. 
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TABLE I: The critical temperatures determined by both Monte Carlo (MC) simulations and also mean-field (MF) 
calculations for Ising films. Bulk 3d results are also given for comparison. Bulk* refers to the rc(oo) parameter of 
the fitting expression of eqn. (jH)). Critical exponents are given which have been extracted from the MC simulations 
via a finite size scaling analysis. Since we can extract the ratio of exponents ^/v for the bulk 3d system we give 
(7/1/)' as 7/1/ + 2 - 3 
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TABLE II: Fitting parameters for Ising thin- films using eqn. ©. MC and MF stand for Monte Carlo and 
mean-field respectively. 
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a) total magnetisation ^ b) total susceptibility 




monolayer 8 layers a monolayer 8 layers a 

2 layers x ^ 20 layers ^ ■ ^ 2 layers x ^ 20 layers ^ ■ ^ 

4 layers « bulk - o 4 layers * bulk - o 



c) layer magnetisation d) layer susceptibility 




123456 123456 
koT/J koT/J 



FIG. 1: The magnetic properties of Ising sc thin-films from Monte Carlo (MC) simulations, (a) and (b) show a 
dimensional crossover of m and x respectively from 2d- to bulk-like for 128 x 128 x / spins (l is the number of 
layers in the film and I = 128 represents the bulk system), (c) and (d) show the layer-dependence of and 
Xk, where A; is a layer index, in 128 x 128 x 10 sc films. Owing to the isotropic exchange interaction and spatial 
symmetry, layers k and 10 — A; have the same properties. Lines are added as a viewing aid. The error bars are 
smaller than point size. 
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FIG. 2: The magnetic properties of Ising sc thin-films from mean-field calculations. As in Fig. ^ (a) and (b) show 
a crossover of m and x from the 2d to the bulk limit {ksTc/J varies from 4 to 6). (c) and (d) show examples of 
the layer dependence of rrik and Xk for 10-layered sc films. 
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FIG. 3: The critical temperatures Tc as a function of thickness I extracted from (a) Monte Carlo simulations and 
(b) mean-field calculations. For mean-field, the values of ksTc/J are 4 for 2d and 6, 8 and 12 for bulk sc, bcc, 
and fee respectively. Lines are added as a viewing aid. The error bars in (a) are smaller than point size. 
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FIG. 4: An example of the extraction of critical ex- 
ponents /3,(7/i/)', and \/v from the slope of a least- 
square-fit (see text) for 10-layered sc films. The 
apparent linear relation supports eqn. @ 
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FIG. 5: Results of critical exponents for all sc, bcc, 
and fee films. In all structures the exponents are 
very close to 2d values. Especially, for thin-films 
{I < 8 approximately), the exponents are consistent 
with 2d results and suggest that thin-films belong 
to the 2d class. The summation {'y/vy + 2(3/v) = 2 
supports our single correlation length ^ assump- 
tion. 
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FIG. 6: The scaling functions for the magnetisation (a) and susceptibihty (b) in 10 layered sc systems with L 
from 100 to 200 calculated by using those critical exponents extracted for 10-layered sc films with L ranging from 
64 to 128. Analogous plots using 2d critical exponents are shown in (c) and (d). In (a) and (b), we find a good 
collapse of the data for all systems with L > 60 including L = 150, 200 which are bigger than those used in the 
exponent extraction. Note, however, that for L close to the thickness / = 10, the collapsing is not found since the 
assumption underpinning our finite size scaling relation (eqn. ^ rests on the condition L ^ I. In (c) and (d) we 
find a slightly worse collapsing than that in (a) and (b) - compare, for example, the collapsing around x = 1 in 
(c) and (a). This confirms the weak I dependence of critical exponents and the deviation from the 2d values when 
the films thicken. 
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FIG. 7: Comparison of the critical expo- 
nent /? extracted from experimental data for 
Ni(lll)/W(110) in [i] and from Monte Carlo sim- 
ulations of Ising thin- films. 
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FIG. 8: Double logarithmic plot of m against t, the 
reduced temperature, for 4-layered Ising sc films. 
The finite size effect is evident once ^ becomes com- 
parable to the size L of the layers of the films. Note 
the data is reliable down to logt ^ —- log L 
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FIG. 9: The effective critical exponent /3eff extracted from 100 x 100 x I sc Ising films where I is the number of 
layers. For / > 4, it is found that around logt = —1 to —1.5, the /?eff curves reach peaks and this indicates how a 
step-like function of /3 of the form described in [3] i.e. Fig. [3 is possible. 
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